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DETERMINATION  OF  NOISE  FIELD  DIRECTIONALITY  DIRECTLY  FROM 
SPATIAL  CORRELATION  FOR  LINEAR,  PLANAR,  AND  VOLUMETRIC  ARRAYS 

INTRODUCTION 

When  an  array  is  located  in  a  homogeneous  stationary  noise 
field,  measurement  of  the  crosscorrelations  between  all  pairs  of 
separated  elements,  at  each  temporal-frequency  of  interest,  is 
the  most  general  second-order  statistical  information  that  can  be 
extracted.  These  spatial  correlations  depend  upon  the  direction¬ 
ality  of  the  surrounding  noise  field,  which  is  the  primary  quan¬ 
tity  of  interest  here.  Instead  of  beamforming  the  element 
outputs,  for  example,  and  trying  to  suppress  the  inherent  side- 
lobes  by  proper  weighting  procedures,  we  want  to  avoid  any  pre¬ 
conceived  notions  about  data  processing  and  go  directly  from  the 
spatial  correlation  to  the  noise  field  directionality  in  as 
direct  and  simple  a  manner  as  possible. 

However,  because  the  noise  field  directionality  is  a  two- 
dimensional  function  of  polar  and  azimuthal  angles,  some  inherent 
loss  or  condensation  of  information  takes  place  with  a  linear 
array  and,  to  a  much  lesser  extent,  with  a  planar  array.  Never¬ 
theless,  we  want  to  preserve  and  extract  the  maximum  amount  of 
information  about  the  noise  field  directionality,  consistent  with 
the  dimensionality  of  the  array  employed,  and  to  minimize  the 
amount  of  data  processing  required. 
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We  begin  by  assuming  the  array  to  be  an  infinite  continuous 
line  in  the  one-dimensional  case,  and  solve  the  integral  equation 
for  the  integrated  (or  collapsed)  noise  field  directionality,  at 
each  temporal-frequency,  in  terms  of  the  spatial  correlation 
along  the  line.  Then,  we  discretize  the  line,  so  as  to  be  an 
equi-spaced  array,  and  determine  the  effect  that  this  limitation 
has  upon  the  estimated  directionality.  Finally,  we  investigate 
the  smoothing  that  is  caused  by  the  practical  requirement  that 
any  physical  array  must  have  finite  length.  Thus,  the  facts  that 
the  spatial  correlation  will  never  be  available  on  a  continuum, 
nor  for  infinite  separations,  are  included  in  the  analysis. 

A  similar  procedure  is  pursued  for  the  two-dimensional  case, 
where  the  planar  array  is  presumed  to  have  equal  spacings  A^  and 
Ay  in  the  x  and  y  dimensions,  respectively.  Again,  the  aliasing 
effects  are  considered,  as  well  as  the  limitation  of  having  to 
employ  a  finite-size  planar  array.  Finally,  in  the  three-dimen¬ 
sional  case,  where  the  problem  is  overdetermined,  a  plausible  and 
efficient  procedure  -cor  collapsing  the  surplus  information  is 
presented,  although  it  is  recognized  that  an  unlimited  number  of 
alternatives  exist. 

Although  it  was  stated  that  the  noise  field  directionality  is 
of  interest,  this  does  not  preclude  the  presence  of  plane-wave 
arrivals,  that  is,  additive  signals  or  interferences  in  the  back¬ 
ground.  In  fact,  the  examples  are  specifically  of  that  type,  for 
these  can  be  considered  as  the  fundamental  building  blocks  of  a 
general  noise  field. 
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Some  related  results  on  this  problem  of  restoring  the  noise 
field  directionality  from  the  spatial  correlation  are  given  in 
11,2,3,4],  but  limited  to  the  line  array.  Specifically,  [1]  gave 
a  least  squares  approach,  starting  from  a  discrete  finite-length 
array.  However,  ill-conditioning  of  the  simultaneous  linear 
equations  for  the  noise  field  directionality  precluded  its  use 
for  more  than  approximately  ten  elements.  This  ill-conditioning 
is  circumvented  here  by  deferring  the  discretization  until  after 
the  integral  equation  is  solved;  this  procedure  for  the  line 
array  was  first  given  in  [4]. 
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CHARACTERIZATION  OF  NOISE  FIELD 

Let  N^(9,4>)  be  the  intensity  of  the  homogeneous  stationary 
noise  field  at  temporal-frequency  f,  arriving  from  direction  e,*>, 
where  0  <  9  <  n,  -n  <  <*  <  n;  see  figure  1.  The  amount  of  power 
received  in  solid  angle  d9  d<*  sine  about  9,$  is 

d9  d<*  sin9  Nf  (  e,  <*)  .  ( 1 ) 

We  call  Nf(9,<*)  the  noise  field  directionality;  the  product 
sin9  Nf ( 0, ♦)  could  be  called  the  plane-wave  density. 


Consider  general  field  point  x^,y^,z^.  Then  if  the  time  of 
arrival  at  the  origin,  of  the  component  from  direction  9,<*,  is 
zero,  then  the  time  of  arrival  at  x^,y^,z^  is 

r,  -  —  —  (x.  sin9  cos*  +  y,  sine  sin*  +  z.  cos9) ,  (2) 

1  c  I  1  I 

where  c  is  the  speed  of  propagation.  Therefore,  the  transfer 
function  at  x^,y^,z^  applied  to  the  arrival  from  direction  e,<*  is 
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H^ff)  -  exp( -i2nf )  = 

\ 

»  exp  i^  (x^  sin©  cos$  +  y^  sinQ  sin<*  +  z 1  cos9)  ,  I 

✓ 

where  wavelength  X  »  c/f. 

The  elemental  contribution  to  the  crosscorrelation  between 
this  arrival  at  x^,y^,z^  and  X2'Y2'Z2'  at  temPoral-f requency  f, 
is  then 


de  d<i>  Sine  Nf{e,4>)  H  x  (  f )  H2(f)  =  de  d  4>  sine  x 


x  Nf  ( 6,  ♦)  exp  i—r  (x  sine  cos<l>  +  y  sine  sin<*  +  z  cos9) 


where  separations 


X1  -  x2 ' 


*1  -  y2 


z  ®  z  2  z  2  • 


If  the  arrivals  from  different  directions  are  uncorrelated,  the 
spatial  correlation  (at  frequency  f)  between  two  points  separated 
by  x,y,z  is  then  given  by  integrating  over  all  angular  space,. 


n  n 


Gf ( x , y , z ) 


de  d$  sine  Nf ( 6, ♦)  x 


0  -  n 


exp  i  — t-  (x  sin9  cos*  +  y  sine  sin*>  +  z  cos0)  . 


The  problem  of  interest  is:  given  spatial  correlation 
Gj(x,y,z)  versus  x,y,z  (or  restricted  slices  of  Gf(x,y,z)),  solve 
for  noise  field  directionality  N^(e,$)  (or  smoothed  versions  of 
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Nf  ( 0,  «>)  )  .  That  is,  invert  integral  equation  (6)  for  noise  field 
directionality  Nf(9,40  or  for  whatever  can  be  determined.  There 
are  three  cases  that  must  be  distinguished,  namely,  linear, 
planar,  and  volumetric  arrays. 
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LINEAR  ARRAY 


It  is  most  convenient  mathematically  to  locate  the  line  array 
along  the  z  axis,  that  is,  x  -  y  -  0 .  Then  the  exponential  in 


independent 

of  <j>, 

.  and  (6)  reduces  to* 

r 

t  r 

i 

f  \ 

G^  (  z  )  » 

de 

d 4»  sine  Nf  (  0,  ♦  )  exp 

i^  z  cos© 

( 

)  - 

-n 

•.  ✓ 

whe  re 


n 


- 

% 

d0  sin©  N f ( 0)  exp 

/  \ 
.  2n 

i— ^  z  cose 

V 

<  J 

( 

Nf  (  0 ) 


n 

d*  Nf(0,4>)  for  0  <  0  <  it 
-n 


(7) 


(8) 


is  the  integrated  or  averaged  noise  field  directionality,  and 
Gg(z)  is  the  one-dimensional  spatial  correlation  at  separation  z 
along  the  line,  both  functions  evaluated  at  frequency  f.  G^(z) 
is  the  only  second-order  function  that  can  be  measured  (or 
estimated)  from  the  line  array,  and  Nf(0)  is  the  only  field 
runction  that  can  be  determined.  There  is  no  possibility  of 
undoing  the  integration  of  (8);  this  is  a  mathematical  represen¬ 
tation  of  the  inherent  conical  symmetry  of  response  of  a  linear 


*The  case  where  the  line  array  is  located  on  the  x  axis  is 
treated  in  appendix  A. 
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array.  It  is  also  one  reason  for  choosing  the  line  array  to  lie 
along  the  9*0  axis,  since  all  the  two-dimensional  field  infor¬ 
mation  is  conveniently  _ollapsed  into  a  one-dimensional  function 
of  0  alone.  See  appendix  A  for  the  problems  associated  with 
choosing  a  different  coordinate  system. 


SOLUTION  OF  INTEGRAL  EQUATION 

To  solve  integral  equation  (7)  for  noise  field  directionality 
N^(0),  consider  the  following: 

+® 


J  dz  exp^-i^  u  zj  Gf(z) 


+  eo 


*  J  de  sin0  Nf(0)  J  dz  exp^-i^y  z  (u  -  cos©)j  * 
0  -® 


d0  sin©  N^{0)  X  S(u  -  cos©) , 


where  6  is  the  delta  function.  Now  let  t  *  cos©,  which  is  a 
one-to-one  transformation  for  0  <  9  <  n,  to  get 
+<*>  1 

J  dz  exp^-i^y  u  zj  G^(z)  -  X  J  dt  N^(acos(t))  5(u  -  t)  = 
-00  -1 

f  X  Nj(acos(u)  )  for  |u|  <  1  'j 

I  0  for  | u  j  >  1  / 

That  is, 


(9) 


( 10) 
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+<*» 

Nf(acos(u))  =  ^  J  dz  exp^-i^  u  zj  Gf(z)  for  |u|  <  1,  (11a) 

—  00 

or 

+  a> 

Nf  ( 0 )  -  ^  J  dz  exp^-i-^y  cose  zj  G^(z)  for  0  <  0  <  n.  (lib) 

—  ao 

Here,  acos  is  the  principal  value  inverse  cosine  function. 

Compare  (lib)  with  starting  point  (7). 

Thus,  given  the  spatial  correlation  G^(z)  for  all  possible 
separations  z  along  the  line  array,  the  integrated  noise  field 
directionality  is  available  via  a  single  one-dimensional 
Fourier  transform. 
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EXAMPLE 

An  example  is  informative  at  this  point.  Let 
Nf  (9)  -  &(  9  -  9  )  ,  0  <  9q  <  n. 

Then  the  spatial  correlation  is,  from  (7), 

G^(z)  «  sinG 

Observe  that  as  9  ->  0  or  n,  that  is,  endfire  of  the  line  array, 
the  strength  of  this  quantity  decays  to  zero,  due  to  the  sine 
term  in  the  area  element  in  (1).  Substitution  of  correlation 
Gf(z)  into  (lib)  yields  noise  field  directionality 

Nf(9)  -  sin9Q  S(cos9  -  cos0Q)  for  0  <  9  <  n. 

Now  the  delta  function  here  is  located  at  9  »  9  and  has  area 

o 

l/sin0Q.  Thus,  N^( 9)  is  S(9  -  9q ) ,  as  it  should  be;  however,  the 

trigonometric  form  shows  N^(9)  as  the  product  of  two  terms,  the 

first  of  which  tends  to  zero  as  9  -*  0  or  n,  and  the  second  of 

which  has  an  area  that  tends  to  infinity  as  9  -»  0  or  it.  This 

o 

behavior  will  re-occur  in  the  following  investigations. 

We  have  employed  the  following  useful  property  above:  if 
g(x)  has  an  isolated  zero  at  xQ,  then  in  the  neighborhood  of  xQ, 

«(9U))  .  s(g'<x0>  (*-x0>)  -  j g,(x  )|  MX  -  x0)  . 

That  is,  the  area  of  the  delta  function  at  xQ  is  equal  to  the 
reciprocal  absolute  slope  of  the  argument  at  xQ,  if  nonzero. 


exp 


Z  C0S9 


,)• 
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DISCRETE  INFINITE-LENGTH  ARRAY 


If  samples  of  spatial  correlation  G^(z)  at  increment  A  in  z 
are  available,  an  approximation  to  (11a)  is  afforded,  for 
( u |  <  1,  by 

+« 

Nf(acos(u))  =  ^  exp^-i^-  u  A  nj  (  A  n),  (12) 

n»-® 

the  right-hand  side  of  which  has  period  X/A  in  u.  Since  the 
integrated  noise  field  directionality  in  (11a)  is  defined  on  an 
interval  of  length  2,  that'is,  -1  <  u  <  1,  aliasing  will  occur  in 
approximation  (12)  unless  A  <  X/2.  Thus,  the  spacing  A,  between 
samples  of  G^(z),  must  be  less  than  a  half-wavelength  at  the 
temporal-frequency  f  of  interest.  This  is  presumed  true  hence¬ 
forth. 

Now  if  u  is  restricted  to  the  values 


u 


m 


m  X 
M  A 


for  -  |  <  m  <  |  -  1, 


(13) 


which  cover  a  full  period,  there  follows,  for 


m  X 
M  A 


<  1/ 


+  0D 

exp( -i2  Jimn/M )  G^(A  n).  (14a) 

n--® 

The  sum  on  the  right-hand  side  can  be  accomplished  via  an  M-point 
fast  Fourier  transform  when  collapsing  is  employed  [5;  p.5].  The 
resultant  angles  at  which  N^(9)  is  available  are 


N, 


es))  = 


6  -  acos 

m 


f-  -1 

Im  aJ  ' 


or  cos9  -  j  for  -  j  < 
m  n  A  2  - 


s  M  , 

m  <  ■=■  -  1 , 


(14b) 
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provided  that  ~  <  1.  These  values  are  equally  spaced  in  cos0 
space . 

The  right-hand  side  of  (12)  can  be  rewritten  in  the  form 
[5;  pp.  3-4] 


+® 


i  J  dz  exp (-1*1  u  z) 


Gf  (z) 


+  ® 

L 

n=*-“ 


S  (  z  -  A  n ) 


+  «  +  a> 

=  Nf(acos(u))  ©  ^  S  ^u  -  ^  |acos  ^u  -  j  ,  (15a) 

n=»-®  n=-® 

where  ©  denotes  convolution.  The  separation  of  these  aliased 
lobes  (for  n  *  0 )  is  A/A  on  the  u  scale;  then,  since  the  extent 
of  Nf(acos(u))  is  2  on  the  u  scale,  overlapped  aliasing  lobes  do 
not  occur  if  A  <  A/2.  This  is  a  mathematical  back-up  to  the 
claim  under  ( 12 ) . 
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DISCRETE  FINITE-LENGTH  ARRAY 

The  effect  of  a  finite-length  array  can  easily  be  incorpora¬ 
ted  by  modifying  (15a),  so  as  to  include  weighting  w(z).  Then, 
we  have,  for  the  estimated  noise  field  directionality. 


Thus,  not  only  is  the  noise  field  directionality  aliased  at  sep¬ 
arations  A/A  in  u,  but,  in  addition,  it  is  smoothed  by  window  W. 
Sampling,  per  se,  does  not  distort  the  estimated  directionality, 
if  done  finely  enough,  that  is,  A  <  X/2 .  However,  the  finite 
length  of  the  array  always  causes  smearing,  with  a  window  width 
of  the  order  of  A/L  ,  where  L  is  the  effective  length  of 
weighting  w( z ) . 


15/16 

Reverse  Blank 
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is  the  sum  of  the  elemental  components  in  symmetrically- 
arriving  rays  on  opposite  sides  of  the  planar  array;  recall  that 
©  -  n/2  now  corresponds  to  the  plane  of  the  array.  Spatial 
correlation  G^(x,y)  is  the  only  function  that  can  be  measured  (or 
estimated)  from  the  planar  array,  and  Nf(9,<H  is  the  only  field 
directionality  function  that  can  be  determined.  There  is  no 
possibility  of  undoing  the  summation  of  (17);  this  is  a  mathemat¬ 
ical  representation  of  the  inherent  two-sided  symmetric  response 
of  a  planar  array.  It  is  also  one  reason  for  choosing  the  planar 
array  to  lie  along  the  9  -  n/2  plane,  since  the  totality  of  the 
two-dimensional  field  information  is  conveniently  collapsed  into 
a  one-sided  function  of  9,  that  is,  0  <  9  <  n/2. 
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SOLUTION  OF  INTEGRAL  EQUATION 

Consider  the  two-dimensional  Fourier  transform  of  (16), 

+  00 

I ( u , v )  S  ||  dx  dy  exp^-i^y  (ux  +  vy)j  Gf(x,y)  » 

-.CD 

(18) 

n/2  n 

=■  d8  |  d<|>  sin©  ( 0 ,  ♦ )  X2  6(u  -  sine  cos$)  S(v  -  sin0  sin«j>). 

-ii 

Let 

a  -  sine  cos<f>,  6  -  sine  sin$  for  0  <  0  <  n/2,  ~n  <  <J>  <  n.  (19) 


These  relations  can  be  inverted  by  using 


to  give 

sine 


a  +  i0  *  sine  exp(i$), 


|a  +  iS|  -  (a2 


♦  =  arg( a  +  ip)  . 


(20) 


(21) 


Thus,  (19)  is  a  one-to-one  two-dimensional  transformation  in  the 
ranges  0  <  9  <  n/2,  -n  <  $  <  n  allowed  in  (18).  From  (19)  and 
(21),  the  Jacobian  is 


3(a,S) 

cose 

cost 

-sine 

sin$ 

3  ( e, ♦ ) 

cose 

sin$ 

sin0 

cos$ 

-  sine  cose 

-  (■ 

=2  ♦  62)15 

2 

a  - 

Substitution  of  these  results  in  (18)  yields 


(22) 
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-H 


I{u,v)  -  X2  JJ  da  d0  F(a,f3)  (l  -  a2  -  02]  S(u-a)  S(v- 


8) 


i  X2  F(u,v)  (l  -  u2  -  v2] 


r  2  2  x  ' 

for  u  +  v  <1 


othe  rwi se 


(23) 


where  is  a  circle  of  radius  1  located  at  the  origin,  and 

F ( a, 0 )  ■  N^(asin(|a  +  i0|),  arg(a  +  i0)).  (24) 

Here,  asin  is  the  principal  value  inverse  sine  function.  From 
(23),  .(  24),  and  (18),  the  noise  field  directionality  is 


Nf(asin(|u  +  iv|),  arg(u  +  iv))  » 


fr  -  °2  -  -2) 


+  09 


JJ  dx  dy  exp^-i^y  (ux  +  vy)j  Gf(x,y) 


for  u2  +  v2  <  1 .  ( 25 ) 


An  alternative  form  is  available  by  letting 


u  -  sin0  cos$,  v  -  sin©  sin$  for  0  <  0  <  ix/2 ,  -n  <  4>  <  it,  (26) 


namely 


+  09 


Nf  (©,♦) 


COS0 


JJ  dx  dy  exp^-i^y  sin©  (x  cos<f>  +  y  s  i  n  )  J  Gf(x,y) 

-  CD 

for  0  <  0  <  it/2,  -it  <  $  <  n.  (27) 


It  is  interesting  to  compare  this  form  with  starting  result  (16). 
An  alternative,  when  the  planar  array  lies  in  the  y  -  0  plane,  is 
given  in  appendix  B. 
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BEHAVIOR  NEAR  PLANE  OF  ARRAY 

At  first  sight,  the  presence  of  the  cos©  terra  in  (27)  would 
appear  to  be  a  problem  for  9  =  n/2,  which  is  the  plane  of  the 
array.  However,  the  following  example  illustrates  what  is 
happening;  let 

N*(9,<fr)  -  6(9  -  9  )  6 ( $  -  *  )  for  0  <  9  <  n/2,  -it  <  *  <  it. 

1  o  o  o  o  ~ 

Then  (16)  yields  spatial  correlation 

G,(x,v)  «  sin9  expfi-^-r  sin9  (x  cos$  +  y  sin$  )]. 

The  strength  of  this  quantity  tends  to  zero  as  9  ->  0.  Substitu¬ 
tion  of  this  Gj(x,y)  in  (27)  yields  noise  field  directionality 

N-(9,$)  *  sin9  cos9  6(sin9  cos+  -  sin9  cos<J>  )  x 
to  o  o 

x  5(sin9  sin<fr  -  sin9  sin<J>  ). 

o  o 

By  use  of  the  property 

5 (ax  +  by)  5(cx  +  dy)  *  ' 

it  may  be  shown  that  Nj(9,$)  is  5(9  -  9q  )  6(<fr  -  <J>Q ) ,  as  expected; 
however,  the  trigonometric  form  shows  N^  as  the  product  of  two 
terms,  the  first  of  which  tends  to  0  as  9  -►  0  or  n/2,  and  the 
second  of  which  has  impulses  with  area  which  tends  to  infinity  as 
9q-»  0  or  n/2.  Thus,  the  sin9Q  and  cos9  terms  are  not  a  problem 
since  they  are  compensated  by  multiplicative  terms;  however,  they 
may  lead  to  inaccuracies  in  numerical  computation. 
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DISCRETE  INFINITE-LENGTH  ARRAY 


The  form  in  (25)  gives  the  noise  field  directionality  sum 

N^,  defined  in  (17),  as  a  double  Fourier  transform  of  the  two- 

dimensional  spatial  correlation  function  G^(x,y).  If  samples  of 

Gf(x,y)  at  increments  A  in  x  and  A  in  y,  respectively,  are 
r  x  y 

2  2 

available,  an  approximation  to  (25)  is  afforded,  for  u  +  v  <  1, 
by 


N^(asin(|u  +  iv|),  arg(u  +  iv)) 


(i 


A  x 

y 


+  co  4.  as 

x  ^  ^  exp^-i^J  (u  Ax  k  +  v  Ay  j)j  Gf(Ax  k,  Ay  j).  (28) 

k»-“  j--«® 


The  summation  on  the  right-hand  side  of  (28)  has  periods  X/Ax  in 

u,  and  X/Ay  in  v.  Since  the  sum  N^  is  defined  within  the  circle 
2  2 

u  +  v  <1,  overlapped  aliasing  lobes  will  occur  in  (28)  unless 

A  <  X/2  and  A  <  X/2 ;  that  is,  the  spacings  between  samples  of 
x  y 

Gf(x,y)  must  be  less  than  a  half-wavelength  at  the  temporal- 
frequency  f  of  interest.  We  presume  this  to  be  true  henceforth. 
Now  if  we  restrict  u  and  v  in  (28)  to  the  values 


u 


m  X 
m  M  A 


for 


n  X  , 

v  «  -  —  for 

n  N  A 

y 


M  .  .  M  . 

“  ^  >  in  \  a  *  jl  / 


N  n  N  - 


(29) 


both  of  which  cover  full  periods  in  u  and  v,  respectively,  there 
follows 
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Nf  ( asin(  |  um  +  ivj),  arg(um  +  ivn)  )  = 


m  n 

4-  ao  4.  cd 


n ' 


A^  Ay  ^  ^  '  exp(-i2;tmk/M-i2nn j/N)  G^tA^  k,  Ay  j),  (30) 

k«-°°  j  --® 


provided  that 


U_  +  1V„ 
m  n 


m  X  .  n  X 

M  A  1  N  A 

x  y 


<  1 


(31) 


The  double  sum  in  (30)  can  be  accomplished  as  an  MxN  two-dimen¬ 
sional  fast  Fourier  transform,  when  collapsing  is  employed 
(5;  p . 5 ] .  The  resultant  angles  at  which  noise  field  directional¬ 
ity  ( 9 , $ )  is  available  are 


0  <  0  -  asin 

-  mn 


m  X  .  n  X 

»  Ax  1  N  Ay 


<  - 
-  •>  • 


or 


.  .  fm  X  .  n  X  'j  . 

■“  <  *mn  ■  ar«  Si^iir  i 

x  y 


_  m  X  .  n  X 

sine _ -  -  -  +  i  —  -7 — 

mn  M  A  N  A 

1  x  y 


[Gsu’*(§m 5  • 


(  32) 


(  33a) 


where 


M  ,  m  <  M  _  . 

2  <  m  <  £  i  > 


N  ,  n  .  N  . 
~2-n-2~^  ’ 


(33b) 


but  remembering  that  (31)  must  remain  true. 
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The  right-hand  side  of  (28)  can  be  re-written  in  the  form 
[5;  pp.  3-4] 


+  09 

JJ  dx  dy  exp^-i^y  ( ux  +  vy)j  Gf(x,y) 


x 


+  00  +00 
xAx  ^S(x-Axk)Ay  ^  S(y  -  Ay  j )  « 
k  =  -®  j  =-® 


Nf(asin(|u  +  iv|),  arg(u  +  iv)) 


•  £  s(“  -  r)  •  £  5(v  -  i;)  •  <34a) 


k— • 


The  separations  of  the  aliased  lobes  (for  (k,j)  t  (0,0))  are  X/Ax 

on  the  u  scale  and  X/Av  on  the  v  scale.  Then,  since  the  extent 

~  2  2 

of  the  noise  field  directionality  N{  is  u  +  v  <1,  overlapped 

aliasing  lobes  do  not  occur  if  A  <  X/2  and  A  <  X/2.  This  is  a 

x  y 

quantitative  restatement  of  the  claims  made  in  the  sequel  to 
(28)  . 
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DISCRETE  FINITE-LENGTH  ARRAY 


The  effect  of  finite  lengths  in  the  x  and  y  directions  can  be 
incorporated  by  modifying  (34a)  so  as  to  include  weighting 
w(x,y).  Then,  we  have,  for  the  estimated  noise  field  direction¬ 
ality, 


(i 


+® 


||  dx  dy  exp^-i^y  (ux  +  vy)j  Gf(x,y) 


+  00  +  0D 

*  \  E  MX  '  k)  \  E  S(y  -  Ay  j )  w( x , y )  = 
k«-®  j  »-® 


®  Nf (asin( |u  +  iv| ) , 


arg(u  +  iv) ) 


E  E«( 


(34b) 


where  window 


+« 


W(u,v)  *  ||  dx  dy  exp^-i^  (ux  +  vy)j  w(x,y)  .  (34c) 


Thus,  not  only  is  the  noise  field  directionality  aliased  at 

separations  X/Ax  in  u  and  X/A^,  in  v,  but,  in  addition,  it  is 

smoothed  by  window  W.  Sampling  alone  does  not  distort  the 

estimated  directionality  if  done  with  Ax  <  X/2  and  A^  <  X/2 ;  see 

(34a).  However,  the  finite  lengths  of  the  array  always  smears, 

with  window  widths  of  the  order  of  X/Lx  in  u  and  X/Ly  in  v,  where 

L  and  L  are  the  effective  lengths  of  weighting  w(x,y)  in  x  and 
x  y 

y,  respectively. 
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VOLUMETRIC  ARRAY 


We  now  have  the  full  version  (6): 

n  ji 

Gf(x,y,z)  »  J  d©  J  d<j>  sin©  N f(©,<f>)  x 
0  -rt 

x  expj\^y  (x  sin©  cos<f>  +  y  sin©  sin4>  +  z  cos©)j.  (35) 

However,  since  noise  field  directionality  N^  is  a  function  of  two 
variables,  while  spatial  correlation  has  three  arguments,  some 
of  the  information  in  G^  is  superfluous  and  must  be  reduced  or 
collapsed  in  some  fashion. 


SOLUTION  OF  INTEGRAL  EQUATION 

We  begin  by  defining  triple  Fourier  transform 


+  «D 


I(u,v,w)  a  JJJ  dx  dy  dz  q(z)  exp^-i^y  ( ux  +  vy  +  wz)j  G^(x,y,z) 


it  it 

■  d©  J  d«f>  sin©  Ng(9,$)  X^  6(u  -  sin©  cos<|>)  S(v  -  sin©  sin<|>)  x 
-n 

x  Q( w  -  cos©) ,  (  36 ) 

where  we  use  a  weighting  q(z)  on  the  z  variable,  and  define 

+  CO 

Q(t)  -  J  dz  exp( -i2ntz/X)  q(z).  (37) 

—  OD 

(If  q(z)  -  1  for  all  z,  then  Q(t)  -  X  S(t).) 
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We  now  break  the  right-hand  side  of  (36)  into  two  parts 
according  to 

it  rt  a/2  a  a  a 


a/2  a 

J  d©  J  d+  -  J  de  J  d+  +  J  de  J  d*  , 

0  -a  0  -a  a/2  -a 


(38) 


and  in  each  region,  we  make  the  change  of  variable  used  in  (19) 
et  seq. ,  namely 


a  =  sine  cos<j>,  (3  “  sine  sin<J>. 


(39) 


Then  <J>  -  arg(a  +  i(3),  while 


6 


asin(|a  +  i0|  )  for  0  <  e  <  a/2' 

a  -  asin(|a  +  i0|)  for  a/2  <  9  <  a 


(40) 


and 


f  (l  -  a2  -  02)  for  o  <  e  <  a/2) 


cose  «  ( 


l  -(l  -  /  -  62)  for  a/2  <  6  <  a  J 


(41) 


Define,  for  future  use. 


(a,  6)  -  (l  -  a2  -  02)  for  a2  +  62  <  1 


Using  these  results  in  (36),  there  follows 


(42) 


I(u,v,w)  -  X2  JJ  do  d0  s( o, 8)  ^  F^(a,0)  6(u-a)  S(v-0)  Q(w-s(a,0)) 
C1 

+  X2  JJ  da  d0  s ( a , 8 ) _1  F2(a,0)  S(u-a)  S(v-6)  Q(w+s(a,0)),  (43) 


where  is  a  circle  of  radius  1  located  at  the  origin,  and 
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F 1  (  a,  0) 
F2(a,0) 


Ng(asin(|a  +  i0|),  arg(a  +  i(3)) 

N^()i  -  asin{|a  +  i0|),  arg(a  +  i0)  L 


for  | a  + 


Evaluating  the  integrals  in  (43),  we  have 

2  -1  r 

I(U,V,W)  -  X  S  (  U  ,  V  )  IF^U/V)  Q(w  -  s  ( u ,  v )  )  + 

1  2  2 

+  F2(u,v)  Q(w  +  s(u,v))|  for  u  +  v  <  1. 


SIMULTANEOUS  EQUATIONS 

If  we  evaluate  the  triple  Fourier  transform  in  (36) 
different  values  of  w,  we  have 

S(U^V'  I  (  u,  v,  )  -  Fj.fUrV)  Q(w1-  S  (  u ,  v )  )  +  F2(u,v)  Q(w2  + 
X 

and 

I(u,v,w2)  -  F^UjV)  Q(w2~  s  ( u ,  v )  )  +  ?2  ( u ,  v )  Q(w2+ 

X 

Also,  if  we  define 

Qn(+)  -  Q(wn  +  s(u,v))  , 

Qn(-)  -  Q(wn  -  s(u,v))  , 

and  denominator 

D  -  Q:(-)  Q2(+)  -  Q:(+)  Q2(-) , 
the  solutions  to  (46)  and  (47)  are 


i0|  <  1. 

(44) 


(45) 


at  two 

s(u, v)  ) 

(46) 

s(u,v) )  . 

(47) 


(48) 

(49) 
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F^u,  v) 
F2(u,v) 


s(u,v) 
X2  D 


s(u,v) 
X2  D 


(q2(  +  )  I(u,v#tf1)  -  Q1(+)  I(u,v,w2)jl  for 


[Q;!^  ( — )  I(u,v,w2)  -  Q2(-)  I(u,v,w1)j 


2  2  - 
u  +  v  <  1 , 


(50) 


provided  that  D  *  0 .  Function  s  is  defined  in  (42). 

There  is  a  great  deal  of  leeway  in  these  solutions.  Namely, 
Q(t)  in  (37)  is  arbitrary,  and  the  values  and  w2  are  arbitrary 
as  well;  the  only  restriction  is  that  D  in  (49)  not  be  zero.  In 
the  special  case  where  weighting  q(z)  in  (36)  is  real  and  even, 
then  Q ( t )  in  (37)  is  also  real  and  even;  we  presume  this  to  be 
the  case  henceforth.  If  we  then  choose  w2  «  -w^ ,  (49)  becomes 

D  =  Q2(w1  -  s(u,v))  -  Q2(w1  +  s(u,v)).  (51) 


If  the  effective  length  of  weighting  q(z)  is  L  ,  a  represent- 

z 

2 

ative  plot  of  Q  (t)  is  displayed  in  figure  2.  For  small  s,  a 

2 

good  location  for  w^  is  at  the  point  where  Q  (t)  has  its  maximum 
slope.  For  larger  s,  a  value  for  w^  near  s  would  guarantee  a 
large  value  for  D  in  (51). 
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ANGULAR  REPRESENTATIONS 

If  we  make  the  substitutions 

u  ■  sin©  cos$,  v  -  sin©  sin<J>  for  0  <  ©  <  n/2,  -n  <  $  <  n 
then  (42),  (44),  and  (51)  yield 

s(sin©  cos<|>,  sin©  sin$)  *  cos© 

F^(sin©  cosK  sin©  sin<)>)  *  N^(©,<^) 

p2(sin©  cos<  ,  sin©  sin<j>)  -  N^(ii  -  ©,♦) 

2  2 
D  -  Q  (\*1  -  cos 0)  -  Q  (w^  +  cos 0) 

Then  (50)  becomes 

Q<wi  “  cos©)  I(+)  -  Q(w.  +  cos©)  I(-) 

Nf<e,*>  -  -  , 

A  Q  (wx  -  cos©)  -  Q  (wx  +  cos©) 

Q<W1  -  cos©)  I(-)  -  Q( w,  +  cos©)  I(+) 

M  _  a  x\  COS©  1 _ 1 _ 

A*  Q/w^  -  cos©)  -  Q  (w^  +  cos©) 

for  0  <  ©  <  n/2 ,  -n  <  <P  <  n, 


for 

0  <  ©  <  n/2,  -rt  < 


where  we  define 


+« 


I(  +  )  -  I(sin©  cos<t>,  sin©  sin$,  +w^ )  -  JJJ  dx  dy  dz  q(z 


expj^-i^  (x  sin©  cos$  +  y  sin©  sin$  +  z  w^)j  Gf(x,y,z) 


,  (52) 

4»  <  n . 
(53) 

(54a) 

(54b) 

)  x 

,  (55) 
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upon  use  of  (36).  We  repeat  that  these  results  for  the  noise 
field  directionality  apply  only  for  Q(t)  real  and  even;  other¬ 
wise,  Q(t)  and  w^  are  arbitrary. 

An  alternative  form  to  (54b)  is  available,  if  desired,  by  the 
substitution  ©'  «  n  -  0,  namely 


cose'  Q<wi  +  cos9' >  I'(-)  ~  ~  cos9')  I'(+) 

X2  Q2(w^  -  cos 0' )  -  Q2(w^  +  cos0' ) 

for  ji/2  <0'<ji,  -rc<<J><n,  (56) 


where 

I'(±)  ■  I(sin0'  cos<j>,  sin9'  sin<f>,  ±w^ )  .  (57) 

The  most  extensive  calculation  required  here  is  that  given  by 
(55);  rewriting  it  differently. 


Ku.v.+w^)  -  JJ  dx  dy  exp^-i^y  (ux  +  vy)j 


+  «D 

X  J  dz  exp^-i^y  q(z)  Gf(x,y,z).  (58) 

The  innermost  integral,  the  Fourier  transform  on  z,  only  needs  to 
be  accomplished  for  the  two  values  ±w^,  whereas  the  outer  inte¬ 
grals  must  be  done  for  ranges  of  u  and  v.  This  is  the  collapsing 
operation  alluded  to  under  (35).  On  the  other  hand,  the  inner 
integral  must  be  repeated  for  every  x,y  value  of  interest; 
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nevertheless,  (58)  is 
Fourier  transform. 

An  example  of  thi 
carried  out  in  append 
taken  with  respect  to 


not  as  difficult  as  a  three-dimensional 

s  procedure  for  the  volumetric  array  is 
ix  C;  it  illustrates  the  care  that  must  be 
the  0  variable  in  (54). 
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SUMMARY 

The  noise  field  directionality  for  the  cases  of  one-,  two-, 
and  three-dimensional  arrays  have  been  solved  for,  explicitly,  in 
terms  of  the  appropriate  spatial  correlation  available  in  each 
case.  In  the  one-dimensional  case,  only  the  polar  directionality 
can  be  determined,  while  in  the  two-dimensional  case,  the  sum  of 
symmetrically  arriving  rays  on  both  sides  of  the  planar  array  can 
be  evaluated.  For  the  three-dimensional  case,  all  ambiguity  can 
be  eliminated,  but  the  overdetermined  nature  of  the  problem 
requires  some  collapsing  of  information  and  leaves  many  options 
to  consider.  For  example,  one  could  let  the  volumetric  array  be 
a  thin-shelled  sphere;  however,  the  resulting  two-dimensional 
integral  equation  for  the  noise  field  directionality  cannot  be 
solved  explicitly.  The  attractive  feature  of  large  stacked 
planar  arrays  is  that  it  permits  the  use  of  Fourier  transforms 
and,  therefore,  an  explicit  expression  for  the  noise  field 
directionality  in  terms  of  the  three-dimensional  spatial  corre¬ 
lation.  Also,  Fourier  transforms  are  efficiently  evaluated  by 
the  use  of  fast  Fourier  transforms. 

In  this  investigation,  we  have  presumed  exact  knowledge  of 
the  spatial  correlation  G^(z)  or  G^(x,y)  or  G^(x,y,z),  depending 
on  the  dimensionality  of  the  array  employed.  In  practice,  Gf 
must  be  estimated  from  measurements  made  from  a  physical  array; 
in  this  case,  maximum  advantage  should  be  taken  of  the  station- 
arity  and  homogeneity  of  the  noise  field.  Thus,  for  a  line  array 
of  equi-spaced  elements,  Gj(nA)  should  be  estimated  from  all  the 
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available  pairs  of  elements  that  have  separation  nA  in  space  and 
over  the  total  available  observation  time  that  data  have  been 
recorded  on  all  elements. 

A  comparison  [6]  is  underway  between  the  methods  of  this 
report  and  the  Fourier  series  method  given  in  [4],  at  least  for 
the  line  array.  Results  are  similar,  but  not  identical;  in 
particular,  the  aliasing  of  the  Fourier  series  method  is  more 
severe  than  for  the  Fourier  integral  approach. 
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APPENDIX  A.  ALTERNATIVE  LOCATION  OF  LINEAR  ARRAY 

If  we  locate  the  line  array  along  the  x  axis,  that  is, 
y  -  z  ■  0,  then  (6)  reduces  to 


Gf  ( x ) 


J  d9  J  d<|>  sin9  Nj(9,$)  exp  ^i-^y  x  sin9  cos<J>j 


-n 


n/2  n 

1  d9  d$  sin9  N^(9,^)  exp^i^y  x  sin9  cos<f>j  ,  (A-l) 


where 


Nf  (  0,  ♦  )  *  Nf  (  9 ,  $  )  +  Nf(Tl  -  9 ,  $ )  +  Nf  (  9,  -♦  )  +  Nf(rt  -  9,-+)  (A-2) 


for  0  <  8  <  n/2,  0  <  $  <  n.  Therefore,  Fourier  transform 


+<*> 


I(u)  *  J  dx  exp^-i^y  u  xj  Gf(x) 


n/2  n 


d9 


d<j>  sin9  N^(9,<j>)  X  6(u 


-  sin9  cos<t>) 


(A-3) 


Now  let 


s  -  sin9,  t  =  cos<f>, 


( A- 4  ) 


which  are  one-to-one  transformations  in  the  ranges  allowed  in 
integral  (A-3).  Then 
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1  1 

I(u)  -  xf  - — - r  - — - r  N,  (  asin(  s  )  ,  acos  ( t )  )  s  S(u-ts). 

0  (l  -  s2)^  -1  (l  -  t2)* 

( A- 5 ) 

The  innermost  integral  on  t  yields 


'  N^(asin(s),  acos(u/s) ) 

(i  -  uW? 

0 


for  | u |  <  s 


for  | u |  >  s  J 


( A- 6  ) 


thereby  giving 


I(u) 


X  F  - ^ - r  Nf(asin(s),  acos(u/s) ) 

,ii  a  -  *2r  u2  -  u2r 

for  | u |  <  1.  (A-7) 


This  integral  equation  for  noise  field  directionality  Nf  is  more 
general  than  Abel's  integral  equation,  because  limit  u  is  also 

J 

involved  in  one  of  the  arguments  of  Wo  hove  been  unable  to 

simplify  (A-7)  and  extract  any  simple  descriptor  of  the  noise 
field  directionality  analogous  to  (8).  Placing  the  linear  array 
along  the  y  axis,  instead,  encounters  the  same  problem. 
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APPENDIX  B.  ALTERNATIVE  LOCATION  OF  PLANAR  ARRAY 


Suppose  the  planar  array  lies  in  the  x,z  plane,  that  is, 


y  *>  0 .  Then  (6)  becomes 


n  n 


Gf ( x, z ) 


J  d9  J  d<f>  sine  Nf(©,t)  exp(i^y(x  sin©  cost  +  z  cos©)) 


0  -n 


n  ix 


d©  J  d<j>  sin©  Nj(9,<|>)  exp^i^y(x  sin©  cost  +  z  cos©)j  ,  (B-l) 


0  0 


where 


Nf(©,4>)  -  Nf(©,t)  +  N  f  (  © ,  —  t )  for  0  <  ©  <  it,  0<*<rt.  (B-2) 


Then 


I(u,v)  a  JJ  dx  dz  exp^-i^y  (ux  +  vz)j  G^(x,z) 


n  n 


J  d©  J  dt  sin©  N^(0,4>)  X2  S(u  -  sin©  cost)  &(v  -  cos©).  (B-3) 


0  0 


Now  let 


s  -  sin©  cost,  t  *  cos©, 


B-4  ) 


for  which  the  Jacobian  is 


IH4I  -  sin2e  sin*  -  (l  -  t2)15  (l  -  s2  -  t2)’1 

o  l  ,  9  ) 


(  B-5 ) 


Then 
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I{u,v)  *  X2  JJ  ds  dt 


Nf  I  acos ( t ) , acos 


(s(l-fc2)  1) 


(l  -  s2  -  t2) 


& ( u- S )  S(v-t) 


'  X2 


(!  -  u2  -  v2) 


N , 


h  -£ 


/ 

f  u  ) 

\  N 

acos ( v ) ,  acos 

1< 

for  u^  +  v  <1 

\ 

K  (l  _  V  2  )  , 

/ 

0  otherwise 


(  B-6 


Thus,  we  have  the  explicit  representation  for  the  noise  field 
di rectionality , 


/ 

f  u  'i  'i 

acos ( v ) ,  acos 

.(l-v2)  Jy 

dz 


X 


x  exp^-i^y  (ux  +  vz)j  G^(x,z)  for  u2  +  v2  <  1  .  (B-7 

If  we  now  let 

u  »  sin©  cos v  =  cosG,  ( B-8 

this  becomes 

+  OB 

N  ^  (  9 ,  <H  ”  2  liM  JJ  d*  dz  exp^-i^y(x  sinG  cos$  +  z  cos9)j  x 


x  Gf ( x , z )  for  0  <  0  <  n,  0<*<n.  (B-9 

This  is  a  viable  alternative  to  (27).  Compare  with  starting 


result  ( B-l ) . 
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APPENDIX  C.  EXAMPLE  FOR  VOLUMETRIC  ARRAY 
Let  the  noise  field  directionality  be  given  by 
Nf(©,*)  -  5(©-0o)  S ( *-*Q ) »  0  <  0Q  <  n,  -n  <  *Q  <  n.  (C-l) 

Notice  that  arrival  angle  0Q  can  range  over  an  interval  of  length 
n.  We  distinguish  two  cases: 

A:  0  <  3Q  <  n/2, 

B:  n/2  <  0^  <  n.  (C-2) 

o 

From  (35),  the  three-dimensional  spatial  correlation  is 

G,(x,y,z)  -  sin©  expfi^Y  (x  sin©  cos*  + 
r  o  \  a  o  o 

+  y  sin©  sin*  +  z  cos0  )|  .  ( C—  3 ) 

The  problem  addressed  here  is  the  reestablishment  of  (C-l)  by 
means  of  the  solution  procedure  given  in  (54)— (57).  Recall  that 
Q ( t )  is  real  and  even. 

First,  substituting  (C-3)  in  (55),  there  follows 

2 

I(+)  -  sin©  X  $(sin©  cos*  -  sin©  cos*  )  x 
—  o  o  o 

x  S(sin©  sin*  -  sin©  sin*  )  Q(+w.  -  cos©  )  .  (C-4) 

o  o  —  I  o 

Now,  when  we  recall  that  0  is  limited  to  ( 0 , n/2 )  in  (54),  the 
delta  functions  in  (C-4)  are  located  at 

A:  ©  -  ©Q  ,  *  ■  *Q  , 

or  B:  ©  -  n  -  ©Q,  *  -  *Q  .  (C-5) 
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By  means  of  the  two-dimensional  transformation  employed  in  (19)- 
(22) ,  we  find  that 

A:  I(±)  -  coie-  S(e  -  V  *<♦  -  *0'  Q(wl  1  coseo>  ' 

o 

X2 

B:  I(±)  “  1  cos  9  '  r  S  ( 0  -  it  +  9o)  6(4>  -  *Q)  Q(w1  +  COS0O)  .  ( C-6 ) 

I  O  I 

Substitution  of  (C-6)  in  the  numerator  of  (54a)  yields 
X2 

A:  cose  cose”  S(e  -  V  S(*  -  V  C(9'9o)  ' 


B: 


COS0 


loose  I 
o 


6(9-  re  +  eQ)  &(  +  -  4>o)  c(efeo)  , 


(  C-7  ) 


where 


C(  0, 0Q )  -  Q(w1-cos0)  Qfw^-cose^  -  Q(  w^  +  cos0)  Q(w1  +  cos9CJ).  (  C-8  ) 
But  since 

2  2 

c<eo,eo>  *  Q  -  cos 0Q )  -  Q  ( w1  +  coseo)  , 

cu-eo,eo)  *  0,  ( c-9 ) 

we  find  that 


Nf  (  0,  ♦  ) 


6(9  -  9  )  6($  -  <j>  )  for  case  A  't 
o  o 

0  for  case  B  , 


(C-10) 


On  the  other  hand,  substitution  of  (C-6)  in  the  numerator  of 
(54b)  yields 
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R:  cosS  5(8  -  V  5<*  -  V  D,9'V  - 

O 

B:  cos©  |-"e~y  S(9  -  n  +  0Q )  6(*  -  *Q)  D(9,90)  ,  (C-ll) 

o 

where 

D(9,9  )  -  Q(w^-cos9)  Q(w^+cos9q)  -  Q(w^+cos9)  Q(w^-cos6q).  (C-12) 
But  since 

D(eo,eo)  -  o, 

D(rt-90,9o)  -  Q2  ( w1  +  cos90 )  -  Q2(w1-coseo)  ,  (C-13) 

we  find  that 

(0  for  case  A  "> 

^  .  ( C-14 ) 

&(0-n+e  )  $(♦-<!>)  for  case  B  J 
o  o 

This  last  case  could  be  written  in  a  form  similar  to  (56)  as 

N-(9',*)  -  S(9'-9  )  5(*-*  )  for  n/2  <  9'  <  n.  (C-15) 
t  o  o 

In  any  event,  (C-10)  and  (C-14)  confirm  starting  result  (C-l)  for 
the  noise  field  directionality. 
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